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This  report  presents  a complete  treatment  of  the  second-order  statistics  of 
intensity  under  conditions  of  strong  scattering.  A three-dimensional  powerlaw 
striation  model  is  used  in  the  analysis.  For  completeness,  the  special  forms 
for  a strictly  two-dimensional  medium  are  also  derived.  Asymptotic  formulas 
for  the  fade  coherence  time  are  derived.  Simple  asymptotic  results  can  be 
derived  only  for  three-dimensional  spectral  indices  less  than  4.  However, 
naturally  occurring  striations  seem  to  fall  in  this  range,  and  the  simple 
formulas  are  generally  applicable.  ^ 
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I EXECUTIVE  SUMMARY 
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In  severely  disturbed  propagation  environments,  the  scintillation 
structure  asymptotically  approaches  a form  that  can  be  characterized  by 
simple  formulas.  For  example,  the  intensity  statistics  are  accurately 
approximated  by  a Rayleigh  distribution.  Since  the  performance  of  both 
coherent  and  non-coherent  systems  in  a Rayleigh  fading  environment  is 
critically  dependent  on  the  fade  coherence  time,  the  second-order 
statistics  must  also  be  specified  to  complete  the  phenomenology-systems 
analysis  link. 

This  report  presents  a complete  treatment  of  the  second-order 
statistics  of  intensity  in  a power-law  scattering  medium.  A three- 
dimensional  striation  model  is  used,  so  that  the  propagation  angle  de- 
pendence in  a three-dimensional  anisotropic  medium  can  be  properly  evalu- 
ated. For  completeness,  however,  the  special  forms  that  apply  to  a 
strictly  two-dimensional  medium  are  also  included.  Thus,  the  results  of 
the  study  are  applicable  to  the  Wideband  satellite  data  base  as  well  as 
the  various  numerical  simulations  that  have  been  performed. 

An  important  finding  from  this  study  is  that  the  asymptotic  behavior 
of  the  second-order  statistics  of  intensity  under  strong  scatter  condi- 
tions is  critically  dependent  on  the  power-law  index.  The  theory  admits 
three-dimensional  striation  spectral  densities  of  the  form  q a where 
2 < a < b.  When  a < 4,  the  contribution  of  large-scale  structures  is 
strongly  suppressed  by  Fresnel  filtering.  When  a S 4,  the  large-scale 
structures  are  only  weakly  suppressed,  and  that  the  strong-scatter  be- 
havior begins  to  show  some  characteristics  of  scattering  in  a medium  that 
contains  a single  dominant  scale  size. 

When  a < ^ > the  second-order  statistics  asymptotically  approach  a 
form  that  is  independent  of  the  Fresnel  radius,  and  a simple  expression 
for  the  fade  coherence  time  can  be  obtained.  In  light  of  data  from  the 
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Wideband  satellite,  it  appears  that  the  effective  three-dimensional 
spectral  index  is  indeed  somewhat  less  than  4.  Thus,  the  simpler  asymp- 
totic results  are  applicable.  When  » 2 4,  the  intensity  statistics  re- 
tain a dependence  on  the  Fresnel  radius,  and  simple  limiting  forms  cannot 
be  obtained. 

A detailed  summary  of  the  results  is  given  in  Section  VI  as  is  a 
comparison  with  earlier  results  based  on  a gaussian  spectral  density 
function.  Under  conditions  of  weak  scatter,  the  results  converge  to  the 
weak-scatter  forms  that  were  independently  derived  in  Rino  and  Matthews 
(1978).  The  intermediate  range  between  weak  and  strong  scattering  must 
be  treated  by  numerical  computations  or  more  elaborate  asymptotic  ap- 
proximations. 
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II  INTRODUCTION 


A full  treatment  of  multiple  scattering  in  an  extended  three- 
dimensional  medium  is  now  well  formulated  [see,  for  example,  Rino  (1978) 
and  the  references  cited  therein].  However,  it  is  impractical  to  solve 
the  necessary  vector  differential  equations  numerically  except  in  certain 
special  cases.  Thus,  simplified  models  and/or  asymptotic  results  are  in- 
variably used  to  interpret  data.  Fante  (1975)  has  reviewed  the  various 
asymptotic  results  that  have  been  applied  to  laser  beam  propagation  where 
isotropic  scattering  is  an  appropriate  idealization. 

Ionospheric  scattering,  by  comparison,  is  complicated  by  the  fact 
that  the  irregularities  are  highly  anisotropic.  Moreover,  the  ionospheric 
outer-scale  cutoff,  q^,  is  sufficiently  small  that  the  condition 
q [Z  « 1,  where  Z = \z  sec  0/(4tt)  is  the  Fresnel  area,  holds  for  all 
observing  geometries  at  the  lowest  frequencies  of  interest  (A.  < 3m). 

This  fact  can  be  used  to  greatly  simplify  the  theory. 

In  Rino  and  Matthews  (1978)  (hereafter  Paper  I),  the  weak-scatter 
theory  was  reformulated  to  show  explicitly  the  ramifications  of  the  very 
large  outer  scale  Si  ^ = 2n/q  . In  particular,  a closed  form  analytic  ex- 
pression for  the  intensity  scintillation  index  S.  was  derived  by  taking 
the  limit  of  the  integral  expression  as  q approaches  zero.  The  result- 
ing expression  is  valid  when  q^AZ  « 1,  which  is  a very  good  approxima- 
tion as  noted  above. 

In  this  paper,  the  results  for  the  second-order  statistics  of  in- 
tensity are  extended  to  accommodate  strong  scattering.  The  analysis  is 
based  on  the  formulation  of  the  gaussian  phase-screen  model  developed  by 
Gochelashvily  and  Shishov  (1971).  The  general  results  are  reviewed  in 
Section  III,  where  it  is  shown  that,  to  calculate  the  second-order  sta- 
tistics of  intensity,  only  the  phase-structure  function  need  be  specified. 
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In  Section  IV  of  this  paper,  the  limiting  forms  of  the  general  re- 
sults are  computed  as  q^  approaches  zero,  as  was  done  in  Paper  I for 
weak  scatter.  The  analysis  is  identical  to  that  presented  by  Rumsey  (1975). 
For  completeness,  however,  the  special  forms  for  a one-dimensional  phase 
screen  are  also  included.  The  scattering  from  a power-law  phase  screen 
is  governed  by  a single  "universal"  strength  parameter,  U,  which  combines 
the  effects  of  perturbation  strength  and  distance  from  the  scattering 
medium.  The  relationship  between  U and  the  intensity  scintillation  in- 
dex is  discussed. 

In  Section  V,  asymptotic  results  are  derived  for  weak  and  strong 

scattering.  The  weak-scatter  results  were,  of  course,  already  deduced 

in  Paper  I.  Under  conditions  of  strong  scattering,  the  asymptotic  forms 

are  critically  dependent  on  the  power-law  spectral  index,  v.  [The  three- 

dimensional  SDF  lias  the  form  C q “ for  q » q 1. 

s o J 

When  v < 1.5,  simple  asymptotic  formulas  can  be  derived  for  the 
autocorrelation  function  of  intensity.  Indeed,  the  scintillation 
index  approaches  unity  from  below  and  (II ')  - 1 = exp  {-2D(y)},  where 
D(y)  is  the  phase  structure  function.  Thie  behavior  is  demonstrated  for 
a special  case  of  one-dimensional  scattering  in  which  an  exact  result 
has  been  obtained. 

When  v s 1.5,  however,  the  intensity  statistics  retain  an  explicit 
dependence  on  the  Fresnel  parameter  Z,  and  S,  can  exceed  unity.  This 
particular  finding  is  evidently  new.  In  effect,  the  behavior  of  the 
scattering  for  the  more  steeply  sloped  spectra  (1.5  ^ v < 2.5)  is  transi- 
tional between  that  for  a power- law  environment  in  which  large-scale 
structures  are  strongly  suppressed  by  Fresnel  filtering  (0.5  < v < 1.5) 
and  a medium  dominated  by  a single  scale  size  where  strong  focusing  can 
occur. 


The  form  of  the  three-dimensional  index  q where  a ~ 2v  + 1 is  largely 
historical.  The  corresponding  one-dimensional  spectral  index  is  a - 2 = 
2v  - 1.  For  scintillation  studies,  it  is  the  one-dimensional  phase 
spectral  index,  p,  that  can  be  measured.  In  terms  of  v,  p = 2v.  Thus, 
if  v = 1.5,  j = 4 and  p = 3. 
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In  Section  VI  the  results  are  summarized  and  compared  to  the  earlier 
results  (e.g.,  Salpeter,  1967)  that  apply  strictly  to  scattering  in  a 
medium  dominated  by  a single  scale  size,  as  opposed  to  a power- law  con- 
tinuum of  scale  sizes.  A simple  formula  for  the  intensity  coherence  time 
is  derived  that  fully  accounts  for  angle  effects  in  a highly  anisotropic 
medium. 


Ill  GENERAL  RESULTS  FOR  A GAUSSIAN  PHASE  SCREEN 


Mercier  (1962)  developed  a mathematically  cumbrous  but  general 
formula  for  the  second-order  statistics  of  the  intensity  fluctuations  of 
a wavefield  as  it  propagates  away  from  a phase-changing  screen.  A more 
compact  formulation  was  developed  by  Gochelashvi 1}’  and  Shishov  (1971). 
Alternative  derivations  of  the  Gochelashvily  and  Shishov  result  have  been 
presented  by  Taylor  (1972)  and  Rumsey  (1975).  Here  we  review  only  the 
principal  results. 

In  the  phase-screen  model,  a phase  perturbation,  6$(p),  is  imparted 
to  a wavefield  at  some  plane,  say  z = z . Diffraction  effects  cause  in- 
tensity fluctuations  (scintillation)  to  develop  as  the  wavefield  propa- 
gates beyone  z = z . The  structure  of  intensity  scintillation  is  charac- 
o 

terized  by  the  correlation  function 


Ej(Ap;z>  = (Up,z)  I(~p ' ,z)) 


(1) 


where  I(o,z)  = |u(p,z)|",  u(c,z)  is  the  complex  wavefield  at  z,  and  the 
angle  brackets  denote  mathematical  expectation. 

The  free-space  propagation  of  R^(Ap;z)  from  z ^ to  z can  be  calcu- 
lated from  the  integral  expression  [the  Gochelashvily  and  Shishov  (1971) 
result]: 


Rj(A| 


• m 


M^(c*,Ap-8 ,0; z ) exp  (-ia-9  k/zjda  d = 


(2) 


where  k = 2~/>  and  ~z  - z-z  . The  quantity  M.  (cv^hc/^  >0;^^ ; z ) is  the 

o 4 o 

fourth-order  coherence  function  of  the  wavefield  in  the  plane  z = z , 

o 

evaluated  in  terms  of  the  variables 


-.(2)  _ lr-(l) 
O'  - 2 L P 

-*(  2)  ->(3) 

- P - P 

, “*(4) .. 

+ P ] 

(3a) 

~*(3)  lr-(l) 
a = -[p 

. -(2) 

-(3) 

— *(  4 ) _ 

(3b) 

+ P 

- P 

- P 

-(4)  _ l--*(  1) 
O'  - — L P 

— •( 2 ) , -(3) 

- P + P 

“*(4) 

- P J 

(3c) 

Two  assumptions  are  used  in  deriving  Eq.  (2).  First,  the  scattering 
of  the  principal  components  of  the  angular  spectrum  must  be  confined  to  a 
narrow  cone  about  the  direction  of  the  reference  wave.  Under  this  assump- 
tion, the  propagation  effects  are  governed  by  the  free-space  form  of  the 
parabolic  wave  equation 


5u(p , z) 

u 


-i 


2k 


2 


U(p,z) 


(4) 


where  §u/3^  is  the  directional  derivative  along  the  propagation  path 
(Rino,  1978).  Second,  M.  cannot  depend  on  the  "centroid"  variable 
a*1  ^ = —7 ~ p ^ + p1"^  + p'^  + pV  . That  is,  is  statistically  homog- 

eneous, or  independent  of  where  it  is  measured.  Note  that  this  neither 
implies  no c requires  that  lower  order  moments  be  statistically  homogeneous 

To  evaluate  M. , the  phase-screen  model  is  used  to  obtain 

M4(p(  * ',pU),s'  J\p(  l);z0)  = v/oxp[i[6cs(p1)  - £0 ( ) + 50  ( P3)  - 6 0(p4)]}>. 

If  60(p)  is  a zero-mean  gaussian  field,  then  Eq . (5)  can  be  evaluated  by 
using  the  well-known  result 

/x^\ 

(exp  {x}>  = exp  {-  } . (6) 

The  final  form  of  Eq.  (5)  can  then  be  conveniently  written  in  terms  of 
the  phase  structure  function, 

D(Ap..)  = <[50(p(i))  - 50(p(j))]2)  , (7) 


(5) 
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(8) 


d s 


M4(c(1),p(2),p(3),p(4);zJ  = exp{-|[D(Ap12)  - D(Apu>  + D(Ap14> 

- D(Ap24)  + D(Ap,A)  + D(Ap„)]} 


34' 


23 


By  making  Che  appropriate  substitutions  from  Eq . (3),  it  follows 
immediately  that  M4  does  not  depend  on  c/ ^ . Moreover, 

M,(aU)>3(J),0;z  ) = exp  {D(a<'‘>))  +D(a(J))  - 7b(S('2)  +a('J)) 

*4  0' 

-W2>-3<3>)}  . 


(9) 


Thus,  the  homogeneity  of  M,  follows  from  the  assumption  that  the  structure 

4 

function  depends  only  on  the  difference  variable  Ap^.-  It  will  be  shown, 
however,  that  because  of  the  symmetry  of  the  variables  in  Eq . (9),  Eq . (2) 
is  well  defined  even  if  l)(Ap)  does  not  exist. 

The  spectral  domain  version  of  Eq.  (2)  takes  a particularly  simple 
form.  By  direct  computation 


= f J Kr(Ap)  exp  {i?  • Ap]  dAp 


// 


MA(-^.5;°,*o>  exp  f‘iK  ’ ^ d? 


(10) 


By  substituting  Eq.  (9)  into  Eq.  (10),  one  obtains 

5j.(K)  =//exP  {-g (?,?£•)}  exp  { -iK  * ?}  d§ 


(11) 


where 


g(?,T))  = D(|)  + D(T)  - |d(§  + Tj)  - |d(|  - T)) 


2"'*  ' ' 2 

which  is  the  principal  result  that  will  be  used  in  later  sections. 

To  investigate  the  general  behavior  of  g(§,T|),  the  spectral  repre- 
sentation. 


(12) 
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D(5) 


L - coslq  • §)]  i,  (q)  dq/(2rr)‘ 


(13) 


can  be  used  where  6 (q)  is  the  two-dimensional  phase  SDF.  If  this  form 

60 

is  substituted  into  Eq.  (12),  after  some  straightforward  manipulations, 
one  obtains  the  equivalent  representation 


S(?D 


■ 


(q)  sin“  (T)  * q/2)  sii-T  (§  • q/2) 


jifL 


(2tt)‘ 


(14) 


From  the  form  of  Eq.  (14),  it  follows  that  if  § (q)  a.  q , then  g(§;T]) 

00 

is  well  defined  as  long  as  2 < a < b.  From  Eq.  (13),  however,  it  is 
clear  that  the  structure  function  itself  is  only  defined  for  2 < a < 4. 

To  carry  the  development  one  step  further,  consider  the  special 

case  in  which  |g(£,T) I « 1 for  all  significant  T]  values.  The  necessary 

conditions  are  stated  in  Section  V.  By  using  the  approximation 

exp  [-g]  ->  1 - g,  substituting  from  Eq.  (14),  and  finally  using  the 
2 

identity  sin  x = (l  - cos  x)/2,  it  follows  from  Eq.  (11)  that 


>r(<)  - 6(<)  + 4$  (O  sin' 

t 00 


which  is  the  well  known  weak-scatter  result. 


?/(2k)] 


(15) 


To  summarize,  in  the  gaussian  phase-screen  model,  the  intensity  SDF 
can  be  computed  by  evaluating  a single  integral,  Eq.  (11).  Only  the 
phase  structure  function  need  be  specified  for  the  computation.  However, 
the  diffraction  effects,  in  conjunction  witli  the  size  distribution  of  the 
irregularities,  are  completely  characterized  by  the  function  g(§,K  zVk) 
”Eq.  (12)  , which  is  well  defined  independently  of  an  outer-scale  cutoff 

-Ct 

for  phase  SDFs  of  the  form  q with  2 < a < 6. 


IV  SCATTER  IN  STRICT  POWER- LAW  ENVIRONMENTS 


In  Section  III,  we  showed  that  the  second-order  statistics  of  in- 
tensity are  completely  characterized  by  the  function  as  defined 

by  Eq.  (12)  in  terms  of  the  phase  structure  function.  If  the  power-law 
model  in  Paper  1 is  used, 


where 


* iq) 


a bC 


.+1/2 


C - r ■ L sec  G C 
p e s 


(lb) 


(17) 


and 


q 


+ B ' / 
x y 


( 18) 


The  coefficients  A',  By,  and  C'  are  defined  in  Paper  I (Eqs.  2 ba,  2b b, 

and  26c).  The  three-dimensional  irregularity  SDF  has  the  same  form  as 

Eq . (lb),  except  that  abC  is  replaced  by  C . The  parameters  a and  b are 

P s 

axial  ratios  along  and  transverse  to  tiie  principal  irregularity  axis.  A 
complete  discussion  of  this  model  is  given  in  Rino  and  Fremouw  (1977). 

The  general  form  of  the  two-dimensional  phase  SDF  corresponding  to 


Eq.  (lb)  is 

R.^(y) 


GC  q y 

P |_o_| 

2nr(v  + 1/2)  I 2 1 


v- 1/2 


Kv-l/2(V)/qo'-1 


(19) 


where  the  geometric  factor  G is  defined  by  Eq.  (10)  in  Paper  1.  Hereafter, 
we  shall  consider  only  isotropic  irregularities  (a  = b = 1).  The  general 
case  is  easily  retrieved  by  reintroducing  G and  interpreting  y appropriately. 
By  using  the  small-argument  formula, 

K _1/;(x)  ~ 1/2  f(v  - 1/2) |x/2|'<V'1/2) 


12 
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it  is  readily  shown  that 


Rc  (0)  = 


C T(v  - 0.5) 
— II 

4irr(v  +0.5)  q~ 


For  completeness,  a one-dimensional  phase  screen  for  which 


(q) 


['l + '1 


% ^ l¥l  2 ^ ‘V'C' 

2 


r'  (0)  = 


< r(^) 

2Arf(p/2)  qj 


will  also  be  considered.  The  one-dimensional  model  is  appropriate  for 
propagation  across  highly  elongated  irregularities.  The  relationship 
between  p and  v is  p = 2v. 

Whenever  the  phase  autocorrelation  function  (ACF)  exists,  the  phase 
structure  function  can  be  written  as 


D(Ap)  = 2[R  (0)  - R (Ap) ] 

6$  o0 


Thus,  from  Eqs.  (19)  and  (20), 


C T(v  - 0.5)  1 - 2 

D(y)  = 2tt  T(v  + 0.5) 


v-1/2 


R i /->(q  y)/r(v  - 0.5) 
v-1/2  o 

2v-l 


- *.  •.  \r,x 


d'  (y) 


^ rl^l  !_ 

-fit  "(p/2) 


v (p-O/2 


K(P-l)/2(V)/r( 

,P-1 


p - u 


As  q approaches  zero,  the  correlation  functions  R (y)  and  R (y)  be- 
° 6®  6 0 

come  singular.  However,  for  a limited  range  of  v and  p values,  D(y)  and 

D ; (y ) remain  well  defined. 

The  limiting  forms  of  D(y)  and  D/(y)  as  q approaches  zero  are  not 
1 Lm  ^ 

immediate,  since  0 D(y)  = 0/0.  Thus,  one  applies  L'Hospital's  rule, 

whereby  the  limit  is  obtained  as  the  limit  of  the  ratio  of  the  separate 
derivatives  of  the  numerator  and  the  denominator.  The  results  are 


lim  . lim  £ 

q - 0 - q - 0 2tt;-(v  + 0.3) 

o o 


V V-l/2 

2H~1  Vj/2(qoy)  y 

(2v  - l)qJ"'“ 


_il  2 , (1.5  - v) 

2-7  f(v  + 0.5)  (2v  - [)  22'J~l 


(p-l)/2 

, q >■ 

..  C#  2\-$-\  K, 

Um  D'  (v)  = Um  E 2 y 

qo  “ ° qo  - ° (l’/2)  (n  - l)q‘’" 

o 


2r'^— 


K,  „ >')  y 

p-d 


; (p/2) ( p - i ) 


It  follow  from  Eqs.  (27)  and  (28)  that  : > specti  mcti< 


are  well  defined  as  q approach? - z<.ro,  as  long  as  0.  > < 

o 

1 < p < J. 


. 3 and 


However,  it  was  shown  in  S<  t i in  f ( - iat  ) ’oil  <■  f i ed 

beyond  v = 1.5.  Similarly,  f>Cr  T ) 1;  well  defined  o'  p 3. 
Hence,  the  most  general  limiting  form  cannot  bi  dedu  i tut in 


Eq.  (.27)  or  (28)  into  Eq.  (12).  Two  methods  of  avoiding  this  dilemma  have 

been  discussed  bv  Rumsey  (1975).  The  most  straightforward  procedure  is 

to  apply  the  finite  q^  form  of  D(y)  or  [b/(y)]  and  use  1.'  Hospital ' s 

— * — * 

r K to  evaluate  the  limiting  form  of  g('.T).  Thus, 

«<?.*>  - 1 -jvvor  1 1 - nrhrsr  [Jlv 

+ dV2/2l''I/2  Kv-1/2(V2> 

- |V3/2T'1/J  Kv-1/2(V3) 


lqoy4/2l'-l/J  Kv.i/2<q0y4>  /q02'-1 


A similar  ex- 


where  y = | = | , y,  - |'j,  y j = 1=  - T | , and  y^  = | § + Ti  | . A similar  ex- 
pression can  be  generated  for  g7('. ”).  To  evaluate  Eq.  (29)  and  the 
corresponding  expression  for  g,('"."),  L'llospital's  rule  must  be  applied 
twice.  The  results  of  this  computation  are  summarized  in  Table  1. 


The  results  in  Table  1 for  gl?,')  agree  with  Rumsey's  (1975)  Eqs. 
(18)  and  (19)  lor  0.5  - - 1.5,  but  not  for  1.5  < v < 2.5.  The  limiting 

forms  of  g(r,  ) lor  v < 1.5  and  v > 1.5  cannot  be  simply  combined.  A 
ter;  has  evidently  been  omitted  in  Rumsey's  computation  for  v > 1.5 
( , > 4).  In  my  ca  , , - 1.5  marks  an  important  transition  in  the  scat- 
tering behavior,  which  will  be  shown  in  Section  IV. 

To  illustrate  the  significance  of  these  results,  let 
. - L/(  2v- 1 ) \ . .u . 


? = C(v) 


in  Eq.  (11)  where 


• 


2tt 

r(o.5 

L K L.  J - v; 

2v-  1 

+ v)(2v  - 1)2 

0.5 

< v < 1 

C(v)  = *< 

C 

2tt 

v = 

1 .5 

C 

1(2.5  - v) 

2n 

c 

r(o.5 

+ v)(2v  - 1) (2v  - 2) 22v” 2 

1.5 

< v < 2 

Then,  with  an  appropriate  substitution  for  q, 


*I(<50C(v)1'<2''-l>)  = 


, z/k  C(v) 


1/ (v-0. 5) 


cos  (t, 


q')  dg' 


where 


h(§,Tp  = 


* 2|el2v-1  ± 2|tT|2v-1  t I?  - tT|2n>-L  - II  +7I2"-1 


|2  log  |||  - |t||"  Log  |ij|  + 1/2 1 § - 7)  | 2 log  ||  - 


+ 1/  2 1 1 + Log  ||  + Tl| 


The  form  of  Eq.  (31)  suggests  the  definition 


U = (z7k)J  0,5  C(v) 


= C Z 
P 


v-0. 5 


F(  v) 


.5  C 


(3 


* 1.5 


= 1.5 
( 


(3 


where 


The  reason  for  introducing  the  new  definitions,  Z and  F(v),  will  become 
clear  shortly. 


It  now  follows  that  the  second-order  statistics  of  intensity  are 
completely  characterized  by  the  integral 


rr 

/•?  (v-0,5)-/) 

Kq;u)  =JJc\p 

-h\',U  q / 

I(q;U)  = JJ  cxpj*h\C,U  ’ q /|cos  ('  • q ) d£  , (36) 

Similarly,  for  the  one-dimensional  phase  screen 

l/(q/;U/)  = /exp< -h , U7 ^ q j cos  (Tq  )d C (37) 

completely  specifies  the  second-order  statistics  of  intensity,  where 

u'  = a/k/p-1)/2  c'(p) 

= c'  Z(p‘1)/2  f'  (p)  , (38) 


- (^) 


(p/2)(p-l)2 


1 < p < 3 


c'(p)  = 


P = 3 


P F(5/2  - , 1 ) 

. (p/2)(p-l)(p-2)2P"2 


3 < p < 5 


F'(P)  = 2(p-1)/2  c'(p)/c' 


The  definition  of  h,(C,*>)  is  readily  deduced  from  Table  1. 


r 


I I 
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The  single  parameter  U determines  the  general  behavior  of  the  second- 
order  intensity  statistics.  In  Section  V asymptotic  results  will  be  de- 
rived for  weak-scatter  conditions  (U  « 1)  and  strong  scatter  conditions 
(U  » 1).  Before  doing  so,  however,  it  is  instructive  to  compare  the  U 
parameter  to  the  intensity  scintillation  index  S 

4 

[S,2 * 4  = R (0 ; z)  - 1] 

L 4 I J 

Paper  I showed  that  for  isotropic  irregularities 


s2  = c zv-°'5 

4 p 


2/tt  (v-0.5)  I 


^/v+0.5 


Thus,  from  Eq.  (33), 


S4  = U 


r(M^) 


Ifr. 


r(^y^)F(v) 


Similarly,  for  the  one-dimensional  phase-screen , 


(41) 


(42) 


9 


= u' 


2r(^) 

' 4 

- (p-i)  r(£^)F/  (p) 


(43) 


2 

For  a fixed  value  of  v (or  p),  S,  and  F admit  the  same  functional  dependence 

14  2 

on  C and  However,  the  functional  dependences  of  ST  and  l’  on  v are 

P -* 

radically  different. 

2 

From  Eq.  (41)  it  follows  that  ST  is  a continuous  monotonical ly  in- 

4 

creasing  function  of  v.  This  happens  because  as  the  slope  of  the  phase 
spectrum  steepens  with  C fixed,  the  spectral  intensity  of  wavenumbers 
corresponding  to  the  Fresnel  radius  rapidly  increases  relative  to  the 
smal ler-scale  structures.  The  parameter  U,  by  comparison,  is  sharply 


4 


I 


% 
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discontinuous  at  v = 1.5  [U7  is  discontinuous  at  p = 3].  This  is  il- 

? 

lustrated  in  Figure  1,  where  S^/U  is  plotted  as  a function  of  v. 

In  Section  III  it  was  shown  that  Fresnel  filtering  suppresses  the 
influence  of  large-scale  structures  in  intensity  scintillation  as  long 
as  v < 2.5.  The  behavioral  transition  at  v = 1.5  can  be  interpreted  as 
a boundary  between  regimes  in  which  large-scale  structures  are  strongly 
suppressed  (v  < 1.5)  and  those  in  which  they  are  weakly  suppressed 
(v  > 1.5).  The  asymptotic  behavior  of  R^(Ap;z)  for  large  U is  very  dif- 
ferent in  these  two  regimes,  as  we  show  in  Section  V. 


V ASYMPTOTIC  RESULTS  FOR  LARGE  AND  SMALL  VALUES  OF  U 


Following  the  work  of  Gochelashvily  and  Shishov  (1971),  a number 
of  authors  have  derived  asymptotic  formulas  for  $ K; z)  and/or  R^(Ap;z) 
by  using  series  expansions  of  q(5;T|)  (e.g.,  Taylor  and  Infosino,  1976  ; 
Buckley,  1971a,  1971b).  Difficulties  can  arise,  however,  if  g(§,T)  is 
defined  in  terms  of  an  implicit  structure  function.  For  example,  moments 
do  not  exist  unless  an  inner-scale  cutoff  is  introduced.  To  do  so,  how- 
ever, does  not  give  physically  meaningful  results,  since  the  basic  inte- 
gral [Eq.  (11)]  is  well  defined  independent  of  either  the  inner-scale 
or  the  outer-scale  cutoff. 

Proper  treatment  of  this  problem  demands  that  series  approximations 
be  applied  to  Eq.  (32)  directly,  as  Rumsey  (1975)  did  for  the  special 
case  of  v = 1 (a  = 3).  First  note  that  h(§,0)  = 0 irrespective  of  v. 
Hence,  for  sufficiently  small  q, 


(44) 


Substituting  Eq . (44)  into  Eq.  (36)  gives  the  so-called  low-frequency 
approximation 


Note  that  the  product  i C<  ) is  independent  of  t . 

P P 

I 

5 

[he  low* frequency  app  . »n  j Lid  or  q < 

This  is  essentially  equivalent  l Lasl  , • - hi  v'  re  I 


5*4}r.- 
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when  v = 4/3  (a  = 11/3).  The  situation  is  more  complicated  when  q » qf. 
For  v jt  1.5  and  ~ » 5 the  starting  point  is 

, ( v-0. 5 ) §2  + 2(v-0.5)(v-1.5)(§  • aj2 

h(i,')  i2|§|  + 2 r-^ 


where  a„  = 1"/|7l|.  Ordering  the  terms  in  Eq.  (46)  and  retaining  only 
the  lowest  order  terms  in  §/T  gives  the  result 


■2ki2v-1 


0.5  < v < 1.5 


h(§,") 


l(2v-l)[r  + (2v-3)(|  • a)2^  "2v_3  1.5  < v < 2.5  . (47) 


For  v < 1.5,  the  high-frequency  approximation  takes  the  particularly 
simple  form 


1HF^’U')  =//exp  -2Ur  L cos  ^ 


It  follows  from  Eq.  (48)  that  I.,_(q;U)  converges  rapidly  to  zero  for 

HF 

q » 1,  Thus,  if  U « 1,  the  bulk  of  the  intensity  SDF  is  contained  in 
the  low-frequency  approximation.  For  l'  » 1,  cf  is  very  small,  and  the 
bulk  of  the  intensity  SDF  (except  the  singularity  at  q =0)  is  contained 
in  the  high-frequency  approximation.  Thus, 


U « 1 


Kq;U) 


(q)  + I^qjU)  u » i 


This  jehavior  is  nicely  illustrated  by  the  one-dimensional  model 
with  p = 2,  wherein  I(q;U)  can  be  evaluated  exactly.  Indeed,  for  p = 2 


h'(§,*)  = 


r 


By  substituting  Eq . (50)  into  Eq . (34),  it  is  readily  shown  that 


,,/2|  I 

y I q I 


i 7 (q  ;u')  = ?(q)  + 2 


exp  f - 2 | ? 1 1 cos  (q? ) dF 


= 5 (q  ) + 


8 sin"  (u'~  q2/2  ) 2q  sin  (qF  i q I ) 

2 + 2 
4 4-  q 4 4-  q 


exp 


r-2u/2|qh 


4 + q‘ 


1 - exp  f- 2U q 1 1 


(51) 
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For  U « 1,  and  q small,  the  first  sin“  term  dominates  Eq . (51)  in  agree- 
ment with  the  low-frequency  approximation.  For  U ' >s  1.  the  third  term 
dominates  and  there  is  no  further  dependence  on  U/.  The  second  term  in 
the  square  brackets  is  significant  only  in  the  intermediate  q range. 

By  using  the  high-frequency  approximation  it  is  easily  shown  that 
for  large  U 


Rj,(§)  exp  (-2|§|2^  ^ } + 1 


(52) 


In  terms  of  the  original  variable 


, 2v- L, 

Kj(y)  - 1 - exp  -2C(v)y  } 


= exp  (-  D(y) } 


(53) 


rsec  Eq.  (27)  . It  follows  that  S,  = R,(0)  - 1 approaches  unity  as  U 
L J 4 1 

becomes  arbitrarily  large.  Moreover,  upon  recalling  Bramely's  (193-t) 
result  for  a gaussian  phase  screen, 


(uu'  ) = ex  - — •>( y ) j 


(54) 


it  follows  from  Eq.  (53)  that  (II')  - 1 = |(uu'  )|  , again  as  l’  becomes 
arbitrarily  large. 


\ , 
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If  u = x + iy , where  x and  y are  uncorrelated,  equal-variance 
gaussian  processes,  then  it  is  well  known  that  I = |up  is  Rayleigh 
distributed  and  (II/)  - 1 = |(uu")|“.  Moreover,  the  phase  of  U is  uni- 
formly distributed.  Now,  while  Eq.  (53)  alone  cannot  guarantee  uniformly 
distributed  phase  (modulo  2tt)  and  Rayleigh  statistics  for  intensity,  ex- 
perimental results  (Fremouw  et  al. , 1978;  Rino,  1978)  suggest  that  the 
Rayleigh  model  is  a good  approximation. 

For  v > 1.5,  the  high-frequency  approximation  [see  Eq.  (47)]  is 
considerably  more  complicated.  In  particular,  retains  an 

explicit  dependence  on  U such  that  numeric  computations  must  be  used. 
Intuitively,  however,  it  is  clear  that  as  v increases  from  1.5  to  2.5, 
the  intensity  statistics  tend  to  behave  more  like  the  medium  is  being 
dominated  by  a single  scale  size,  even  though  strictly  speaking  the  outer 
scale  is  not  important  until  v exceeds  2.5. 

When  a single  scale  size  is  dominant,  strong  focusing  can  occur 

(Pisareva,  1958),  which  causes  S / to  exceed  unity.  The  form  of 

■4 

I1[r(q;U)  for  v > 1.5  does  not  exclude  this  from  happening  in  a power-law 
environment,  although  the  conditions  that  lead  to  S,  > 1 have  generated 
some  controversy  (Rurnsey,  1975;  Taylor  and  Infosino,  1976). 


VI  SUMMARY  AND  DISCUSSION 


In  this  paper  the  power- law  phase-screen  model  has  been  generalized 
to  accommodate  strong  scattering.  The  intensity  SDF  $j.(q)  was  derived  in 
terms  of  a single  integral  expression  I(q;U): 


has  the  same  functional  form  as  the  weak-scatter  formula  for  the  scintil 

2 2 

lation  index,  S.  . 'Recall  that  C = r \ (L  sec  9)  C ; where  C is  the 
L p e s s 

strength-of-turbulence  as  defined  by  Eq.  (7)  in  Paper  I.'  In  general 


whereas  for  U « 1 


where  IT  is  defined  by  Eq . (34)  in  Paper  I 


It  wac  shown  in  Section  IV  tiiat  as  long  as  v does  not  lie  in  the 

2 2 

range  1.4  ^ v < 1.5,  then  S > U,  where  S.  is  evaluated  using  Eq.  (58). 

4 4 

Thus,  the  weak-scatter  condition  U « 1 need  not  require  that  S,  « 1. 
This  is  consistent  with  the  experience  of  most  experimentalists,  namely 
that  the  weak-scatter  S^  formula  accurately  reproduces  their  data,  even 
for  moderately  large  scintillation  levels. 


As  v approaches  1.5,  F(v)  approaches  infinity.  The  v = 1.5  point 
marks  an  important  transition  from  a scintillation  behavior  regime  in 
which  the  contribution  of  structures  larger  than  the  Fresnel  radius  is 


strongly  suppressed  (v  < 1.5)  to  one  in  which  the  contribution  is  only 
weakly  suppressed  (v  > 1.5).  When  v < 1.5,  the  phase-structure  function 
for  example,  is  well  defined  in  the  limit  as  the  outer-scale  wavenumber 
approaches  zero. 


In  Section  V an  asymptotic  analysis  of  I(q,U)  was  performed  for 
large  and  small  q values.  For  U « 1,  the  low-frequency  approximation 


(2v+l)  . 2 

sin 


accounts  for  all  significant  Fourier  components.  As  U increases,  the 
high-frequency  approximation  accounts  for  all  the  significant  Fourier 
components,  except  the  delta  function  at  q = 0.  The  form  of  the  high 
frequency  approximation,  however,  is  very  critically  dependent  on  v. 


When  v < 1.5,  the  high-frequency  approximation  is  independent  of  U. 
Thus,  under  conditions  of  strong  scattering,  the  intensity  statistics 
converge;  ultimately  there  is  no  further  change  with  increasing  Z.  When 
this  happens,  a simple  analytic  form  for  $^(q)  cannot  be  obtained.  None- 
theless, the  intensity  correlation  function  takes  the  particularly  simple 


To  i aluate  the  temporal  autocorrelation  function  of  intensity,  y is 

replaced  by  v :t,  where  v is  defined  by  Eq.  (13)  in  Paper  I.  Thus, 
el t cll 

it  is  defined  to  be  the  time  sepal  at  ion  to  achieve  an  intensity  decor- 


relation of  e 


\M 


The  formula  becomes  meaningless  as  v approaches  1.5.  For  v £ 1.5,  no 
simple  formula  such  as  Eq,  (60),  exists  and  numeric  computations  must 
be  employed.  Indeed,  in  that  case  will  depend  on  Z,  as  it  does  under 
conditions  of  weak  scatter,  in  which  there  is  also  no  simple  formula. 

It  is  instructive  to  compare  these  results  to  the  earlier  results 

2 2 

for  a gaussian  phase  screen  with  a gaussian  ACF  of  the  form  exp  {-Ap  /rQ}. 
For  the  gaussian  ACF,  the  rms  phase  cp  ^ can  be  computed  from  the  expression 

222  2 -3/2 

cp  = r X (L  sec  0)  G (AN  > 2rr  ' r . (63) 

o e e o 


Under  conditions  of  weak,  isotropic  scattering  (k  « 1) 

ro 


S“  = 2cp2  1 - cos2  (2\z/(nr^)) 


If  Xz/r  « 1, 
o 


S7  2ci 


2/2\z 


L ' 1 


which  is  similar  to  Eq.  (58).  If,  on  the  other  hand , x'z/r"  » 1, 

S2  =»  2cp2 
b4  - 


Moreover,  Mercier  (1962)  has  shown  that  the  intensity  statistics  are 
Rician  in  this  limit.  By  comparison,  in  a strict  power-law  environment, 
Xz  cannot  be  increased  indefinitely  without  approaching  unity.  Put 
another  way,  in  a strict  power- law  environment  one  always  encounters 
strong  scattering  with  increasing  distance  from  the  scattering  region 
before  S,  saturates  independent  of  "z. 

Under  conditions  of  strong  scattering,  the  relationship  (II ')  - 1 = 
exp  {-2D(y)}  where 


f 1 - exp  {-y2/r(}) 


I- 


% 
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will  ultimately  apply.  However,  it  is  known  that  when  1 and 


1 


(68) 


where  is  the  fourth-order  moment  of  the  gaussian  spectrum,  the  scintil- 
lation index  S,  exceeds  unity  (Salpeter.  1967).  This  is  the  strong  focusing 
effect.  This  phenomenon  lias  a counterpart  in  a power-law  medium.  However, 

since  neither  c 0 nor  Q exist,  the  strong  focusing  condition  cannot  be 
o o 

simply  characterized  as  with  Eq.  (68). 

The  scintillation  behavior  in  a gaussian  ACF  environment  is  con- 
veniently summarized  in  diagram  form,  as  was  done  initially  by  Cohen  et 
al.  (1907).  Figure  2,  the  usual  representation,  is  reproduced  from 
Singleton  (1970).  Regions  I and  II  correspond  to  the  near  and  far  zones, 
respectively.  Region  III  is  the  "fully  modulated"  or  saturation  region 
where  S,  is  unity  and  the  relationship  (II ')  - 1 = exp  { - 2 D( y ) } holds. 

Region  IV  is  t lie  strong  focusing  region  where  S,  > 1. 

4 

In  a strict  power- law  environment,  by  comparison,  there  is  no 
characteristic  scale  size  akin  to  r^,  and  the  magnitude  of  the  rms  phase 
is  irrelevant  to  the  scattering  behavior.  Contours  of  constant  have 
the  form  U = const.,  or,  from  Eq.  (56), 

\j.n  s 

C Z = const.  (69) 
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Thus,  if  one  were  to  make  a diagram  similar  to  the  Cohen-Singlcton  dia- 
gram for  a power-law  medium,  the  phase  strength-of-turbulence,  C , would 

2 P 

replace  cpo  and  the  Fresnel  area  \z/4rr  would  replace  Xz/r^. 

The  degree  to  which  Region  II  is  present  depends  critically  on  the 
spectral  index  v.  If  v < 1.5,  does  not  significantly  exceed  unity. 

For  v £ 1.5,  however,  values  substantially  larger  than  1 do  occur. 

When  this  happens,  the  contours  are  still  given  by  Eq . (69).  Thus,  the 

scattering  behavior  in  a strict  power-law  environment  can  be  similar 

2 

to  that  predicted  by  the  Cohen-Singleton  diagram  for  Xz/rQ  « 1 . 

In  fact,  if  a finite  outer-scale  is  introduced,  then  the  Cohen- 
Singleton  diagram  for  a power- law  medium  is  similar  to  that  of  a gaussian 
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FIGURE  2 CURVES  OF  CONSTANT  S4  ON  A 

PLOT  OF  Xz/i-2  vs.  ^0-  -TAKEN  FROM 
SINGLETON  (1970) 


medium  with  the  outer  scale  playing  the  role  of  r . Because  the  outer 

o 

scale  is  at  best  indeterminant  in  the  ionosphere,  however,  such  an 
extrapolation  is  not  useful  for  data  interpretation. 
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